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We present the results of an experimental study of the nucleation of superconductivity at the 
surface of a single crystal YBs in a tilted dc magnetic field. A recently developed experimental 
technique allowed us to determine H C 3 at each side of the sample as a function of the angle between 
the dc magnetic field and the surface. Experiment shows that the ratio H C 3,III C 2 ~ 1-28 in the 
direction perpendicular to the surface dc field while according to the theory this ratio should be 
equal to 1. This sharp distinction cannot be ascribed to the surface roughness. 

PACS numbers: 74.25.Nf, 74.25.Op, 74.70.Ad 



I. INTRODUCTION 



Shortly after Saint- James and de Gennes' prediction 
of surface superconductivity [1] many experiments have 
confirmed their basic idea. Resistive and permeability 
measurements at low frequencies @, d, 0, H, Q showed 
that in parallel to the surface dc field the ratio HcijH C 2 
is close to the theoretical value 1.695. The value of 
H C 3 depended on the angle 9 between the dc field and 
the sample surface and decreased as this angle was in- 
creased. The critical surface current, defined as the cur- 
rent at which the sample exhibits a finite resistance, also 
decreased with increasing 0. The ac response of sur- 
face superconducting states (SSS) had characteristic fea- 
tures such as the loss maximum in the intermediate fields 
H C 2 < Hq < H C 3, low frequency dispersion, and nonlin- 
earity at very low excitation levels [6j. Concurrently, a 
number of theoretical models for the calculation of the 
critical current in the SSS were published (see, for exam- 
ple, references in [||). Agreement between the theoret- 
ical models and experimental data was not satisfactory. 
Therefore, attempts were made to develop a more sophis- 
ticated model of the SSS, in particular, a model of the 
surface vortices was elaborated @, B Q ■ In this model 
the surface vortices could move only if the surface current 
exceeds some critical value. Actually this is a variant of 
the critical state model which was applied to the SSS. 
When the dc magnetic field is parallel to the sample sur- 
face, the normal component is zero and the density of the 
surface vortices is zero. It was proposed 0] that due to 
the surface roughness, locally the dc magnetic field had 
a normal component which provided some finite density 
of surface vortices. The experiments @, Q demonstrated 
that the surface critical current dramatically depended 
on the surface properties. Losses in the critical state 
model have a threshold character with respect to the 
amplitude of excitation G : ] . All the above mentioned ex- 
periments were performed on polycrystalline metals or 
alloys. In spite of all these researches an adequate theory 



that could explain all observed peculiarities has not yet 
been proposed. In addition, recent experiments did not 
reveal any vanishing of losses in weak ac fields [l(| [H| 
and showed that the nonlinear resp onse could not be de- 
scribed by perturbation theory [ill ]. Resistive and torque 
measurements @,H,[1] demonstrated that SSS did not ex- 
ist (H c3 /H C 2 = 1) if the dc field was perpendicular to the 
surface as it was predicted by the theory [l| . On the other 
hand, ac measurements [i| demonstrated that in a per- 
pendicular field the transition to a superconducting state 
takes place in fields above H C 2 . It was assumed that there 
was a superconducting network of negligible volume that 
was not detected in any bulk effects [4[. The main goal 
of this paper is to demonstrate that SSS could exist at 
H > H C 2 perpendicular to the surface dc field, and this 
effect could not be ascribed to the surface roughness of 
the sample. 

II. EXPERIMENTAL DETAILS 

We have measured the ac response of a YBg single 
crystal with sizes 10 by 3 by 1 mm which was cut from 
a large crystal. Sample surfaces were mechanically and 
then chemically polished. DC magnetization curves were 
measured by a commercial SQUID magnetometer. The 
ac response at the fundamental frequency and the third 
harmonic signal were taken concurrently using the pickup 
coil method [13]. A block diagram of the experimental 
setup is shown in Fig. [1] The ac field is supplied by 
two identical drive coils connected in series with the load 
resistor. The sample is inserted into one of the pickup 
coils. Empty pickup coils are balanced without any ex- 
ternal circuits. A small unbalanced signal of the empty 
pickup coils was subtracted as well. The electromagnet 
was rotated around the vertical axis, Fig.[TJ Small probe 
coil (not shown in Fig. [I} is employed for ac amplitude 
calibration. 

Fourier analysis of the time-dependent magnetization 
in an ac field h(t) = 2ft,ocos(wi), yields an expression 
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dc field inclination angle 9. 



III. EXPERIMENTAL RESULTS 

YB6 is an isotropic BCS superconductor. It has some 
advantages due to experimentally well define H C 2 as one 
can see in the inset to Fig. [2^. Properties of YB 6 were 

measured with high accuracy in Refs. m m where it 

was found that T c w 7.2 K and the Ginzburg-Landau 
(GL) parameter n ss 3.5. 

Fig. [2^ presents in- and out-of-phase components of 
Xi for T = 4.2 K as a function of the dc field for some 
angles 9 between the wide sample side (KZ-plane) and 
the dc field while it remained parallel to the narrow side 
(XZ-plane) . Fig.[2jD shows the corresponding results for 
another orientation of the dc field when it was parallel 
to the wide sides and had an angle with the narrow sides 
(see inset (b) to Figs. [T]). One can readily see that for 
these two dc field orientations (Gl and G2) x'i(Hq) and 
Xx(Hq) are well distinguished. 



FIG. 1: (Color online) Block diagram of the experimental 
setup. LFG - low frequency generator, LIA - lock-in amplifier, 
1 - electromagnet poles, 2 - sample, 3 - pickup coils, and 4 - 
drive coils. Inset (a) - sketch of the dc field rotation in XZ- 
plane. Inset (b) - sketch of the dc field rotation in YZ-plane. 



of the form m(t) = ho J2 n X«( w ) ex P( — inut). The fre- 
quency and amplitude of the applied ac field were 565 Hz 
and 0.05 Oe, respectively. This frequency has been cho- 
sen by taking into account our experimental constrains 
and results of that the frequency dispersion was not 
important for the managing superconducting transition. 
We assume that in zero dc field at low temperatures the 
losses are negligible and the ac in-phase component, %' l5 
does not depend on frequency and equals to — 1/47T. This 
assumption permits us to get the values of xi an d X3 
in absolute units at finite dc fields. The absence of fre- 
quency dispersion in zero dc field was verified experimen- 
tally. 

Two experimental configurations were chosen for the 
parallelepiped shaped sample. The first one, (Gl), when 
the dc field rotates in XZ-plane (Fig. [TJ inset (a)), and 
the second one, (G2), when the dc field rotates in YZ- 
plane ( Fig. [TJ inset (b)). With the first configuration, 
Gl, dc field remains parallel to the small faces (XZ- 
plane) for any angle 9 while with G2 the field is always 
parallel to the wide faces ( YZ-plane) . Angle 9 is the angle 
between axis Z and direction of the dc magnetic field. In 
all measurements the ac field was parallel to the axis Z. 
The measured ac response for these two configurations 
was different and therefore it permitted us to calculated 
the ac response for each sample face as a function of the 



IV. DISCUSSION 

The raw experimental data (Fig. [5^ and HJd) contain 
the contribution of four sample sides to the ac response. 
In order to separate the response from each side we used 
the following procedure. Let consider the ac response of 
the long rectangular slab in the normal state with a thin 
superconducting sheath (of thickness several coherence 
length) at the surface (inset to Fig. (2b) . The ac field in 
the bulk could be found from the solution of the two- 
dimensional equation 



d 2 h z d 2 h z 
+ 



dx 2 dy 



2i 
S 2 



h z = 0, 



(1) 



where 5 — cj\j2-Kau> is the normal skin depth and a 
is a normal conductivity. We took into account that 
an applied ac field has only one ^-component and ne- 
glected the small demagnetization factor along the Z-axis 
« 0.045. A surface sheath is capable of carrying a sur- 
face current J s . This current screens the inner part of 
the sample and, neglecting the thickness of this sheath, 
we can write the boundary condition for Eq.(fT]) at each 
sample side in the form 



h + Air J si /c, 



(2) 



where h sl is the amplitude of ac field on the boundary 
between normal sample core and the surface sheath at 
i-th sample side, and J S i is surface current at i-th side. 
Indexes i = 1 , 2 correspond to the wide sides while i = 
3, 4 to the narrow sides. The straightforward solution of 
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FIG. 2: (Color online) Field dependence of x'l an d Xi f° r 
different angles 0. Panel (a): xi f° r orientation of the dc field 
that is parallel to the narrow side of the sample and has an 
angle 6 with wide side (rotation in XZ-plane) . Magnetization 
curve is presented in an upper inset. 

Panel (b): Xi f° r orientation of the dc field that is parallel to 
the wide side of the sample and has an angle 9 with narrow 
side (rotation in YZ -plane). The inset shows the sample cross 
section. The darkened area shows (not in scale) the sample 
part in a normal, and light colored part in a superconducting 
state, respectively. 



Eq. (P) is 



h z = — sin(7rn/2) cos(Trnx/2L x ) 

z — ' irn 

n=l,3... 

h xl smh[k n (Ly + y)} h s2 smh{k n (L y - y)] 



sinh[2fc n Lj, 
4 



nn 



smh[2k n Ly] 
sin(7rn/2) cos(7rny/2L y ) (3) 



+ £ 

n=l,3... 

f h x3 smh[q n (L x + x)] h s4 sinh[g ra (L 



smh[2q n L x ] 



smh[2q n L x 



where 

kl = (irn/2L x f - 2i/5 2 , q 2 n = (nn/2L y ) 2 - 2i/5 2 , 

L x and L y are sample sizes in Xl^-plane ( L x > L y ). The 
average ac magnetic field, h z = L X L J^" f^ y h z dxdy, is 



hfZt + hfZ 2 



(4) 



where 



and 



Zl =8 y: 



tanh(fc m L y ) 
n 2 m 2 k m L y 

m— 1,3 y 



Z 2 = 8J2 



m— 1,3 



tanh(g m La;) 
ir 2 m 2 q m L x ' 



It was taken into account that due to the symmetry h sl — 
h s2 and h s3 = h sA . The observed magnetic susceptibility 



Xi = (hz/ho - 1)/4tt. 



(5) 



Using these expressions and experimental data one can 
find the surface currents at each side of the sample. If 
Hq lays in the XZ-plane and the angle between dc field 
and QZ axis is 8 



1 + 4ttxi = Zxhf/ho + Z 2 h s //h 



(6) 



while if Hq lays in the FZ-plane and forms the same 
angle with 0Z axis, then 



1 + 47TX! = Zihf/ho + Z 2 h a z x /h 



(7) 



where xi and Xi are susceptibilities for Gl and G2 con- 
figurations of the dc field, respectively. These two equa- 
tions allow us to find both h s z and h s2 from experi- 
mental data and to calculate the surface currents from 
Eq. ([2]). The bulk conductivity in a normal state a can 
be found from the ac response at temperatures T > T c . 
Under these conditions h sl = h s3 = Hq. Fig. [3] demon- 
strates the frequency dependence of x" at T—7.5 K. This 
experiment was carried out on the setup described in 
Rcf. [14]. Mapping these data by Eq. © we obtain 
a = 8 x 10 16 CGS. This value is in good agreement with 
the result of Ref. [l3| obtained for actually the same sam- 
ple (it was cut from the same bulb). 

Experimental results can be presented in terms of the 
surface current, J s , normalized in such way that the 
complete shielding by the surface current corresponds to 
J s = — 1 while for the sample in a normal state J s = 0. 
Figs.HJi and[3j3 shows J s as a function of the dc magnetic 
field for several angles. Fig. [3k. corresponds to the case 
when the dc field makes some angle with the surface while 
Fig. 2)3 refers to the case when Ho and h o are parallel 
to the surface but there is an angle between them. Wc 
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FIG. 3: (Color online) Frequency dependence of x'i f° r the 
sample in a normal state at T=7.5 K. 



see that in the perpendicular dc field the superconduct- 
ing transition at the surface takes place at Hq > H C 2- 
The observed response at the third harmonic confirms 
this conclusion also. The third harmonic vanishes if the 
dc field becomes smaller than H C 2- The nonlinear re- 
sponse disappears simultaneously with the disappearance 
of the surface current in an increased dc field. Moreover, 
the nonlinear response depends on the orientation the dc 
field, as is shown by Figs. EH and Et>- The orientation 
of the dc and ac magnetic fields here was the same as 
for Figs. and [2b- For third harmonic data we did not 
perform the procedure of separating the response from 
each sample side. However, the first harmonic data show 
that the wide sides provide the largest contribution to 
the ac response. This is why the dependence of | ^3 1 on 
the angle shown in Fig. Ek reflects the dependence of the 
nonlinear response by SSS on the angle between a paral- 
lel to the surface ac and dc fields. When Ho is parallel 
to the surface, the observed onset of the superconducting 
transition does not depend on the angle between the ac 
and dc fields as was expected, Fig. HJd, but the value of 
J s depends on this angle. Fig. [5] demonstrates the an- 
gular dependence of the surface current J s for several dc 
fields. This dependence is nonmonotonic at low dc fields. 
Detailed discussion of the angular dependence of J s will 
be reported in our forthcoming publications. 

Existence of the SSS in the perpendicular dc field for 
Hq > H C 2 can be ascribed to the surface roughness. It 
is well known that the properties of the surface strongly 
affects the ac response of SSS and H c % depends on the 
surface roughness too 0, H, H3]. SSS are localized in the 
surface layer with a thickness of a few coherence lengths. 
For a rough surface the dc field is not parallel to the local 
normal even in the perpendicular to the surface dc field 
and one could expect the manifestation of SSS in fields 
Hq > H C 2- The value of H c ^ could be calculated from 
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FIG. 4: (Color online) Panel (a): Field dependence of the 
surface current in the case when dc field makes angle with the 
surface. Panel (b): Field dependence of the surface current 
when both dc and ac fields are parallel to the surface but 
make some angle with each other. 



the first GL equation for a given magnetic field 

ln(T c /T){-* + I*! 2 *} + (iV + ~a ) 2 * = 0, (8) 

with boundary conditions ^j=?\s — at the surface. Here 
"J/ = A/Ao is the dimensionless order parameter, Ao is 
the order parameter at a given temperature and zero dc 
field, ~a — 2eA/hc, A is a vector potential (curl A 
1 1 the length unit is £0 - the coherence length at T = 0, 
and n is a local outer normal to the sample surface. For 
a really rough surface this is a 3D equation and in order to 
simplify the problem we considered only the 2D version 
where the surface profile depended only on one direction 
(Z-axis). Because the analytical solution is not possible 
for a rough surface profile we used the numerical method. 
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FIG. 5: (Color online) Field dependence of \xs\ for different 
configurations of the dc field. Panel (a): the case when the 
dc field remains parallel to the YZ sample surface but forms 
some angle with respect to the ac field. Panel (b): the case 
when the dc field remains parallel to the XZ sample plane. 



FIG. 7: (Color online) Panel (a): The model profile of the 
rough surface and orientation of dc magnetic field. Panel 
(b): Calculated order parameter (absolute value) for SSS for 
the model surface in the dc field H > H C 2 and 8 = 45°). 
Panel (c): Measured by AFM actual profile of the surface. 
Panel (d): The calculated order parameter near the surface 
for actual surface profile (6 = 45°). £o is unit length for (b) 
and (d) panels. See text. 
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FIG. 6: (Color online) Real and imaginary parts of the surface 
current as a function of the angle between ac and dc fields for 
several dc fields. Both fields are parallel to the surface. 



Eq. ((SJ had been transformed to 

ln(T c /T) j-*(n,m) + |*| 2 *(n, m)| - 

\Jr(n + 1, m)ty(n - l,m) - 2*(n, m) 
dx 2 

*(n + 1, m)^f(n - l,m) - 2\E r (n, m) 

dz 2 + 
i\$(n, m + 1) - ^(n,m- l)]mdzH + 

tt2 2 7 2 T / \ d^(n,m) 
Hnmrdz^in, m) = 7 — — — - 
K ' ' dt 

at the 2D-grid with steps dx and dz. The stationary non- 



trivial solutions of this set of equations was looked for, 
and the maximal Hq, for which such solution could be 
found, was considered as H c3 . Actually, a grid with 10 4 
points was used. Calculations were performed for several 
types of the surface roughness, for a model with sinu- 
soidal surface profile, Fig [7^,, and for the actual surface 
profile that was measured by an atomic force microscope 
(AFM), Fig. Eh. 

The shape of the superconducting nucleus near the 
rough surface is shown for illustration by Fig (TJ) and 
Fig UJi. Fig [7)} demonstrates the results of the numeric 
calculation of the spatial distribution of the order param- 
eter near the surface with a sinusoidal roughness with pe- 
riod m x £0 and amplitude n x £0 while Fig. [Jji shows the 
order parameter for the surface profile measured by the 
AFM. The angle between the dc field direction and the 
sample surface equals 45° in both cases. Careful analy- 
sis of these pictures shows that a superconducting phase 
nucleates near such points where the angle between dc 
field and local normal is maximal, close to 7r/2. 

Fig. [5] shows both the experimental and calculated 
H c z/H C 2 (9) dependence for several models of surface 
roughness. Curve (A) corresponds to the ideally smooth 
surface; (B) to a sinusoidal model with period 20 x £0 and 
amplitude lx(o; (C) - for period 10 x £0 and amplitude 
5x( ; and curve (D) correspond to the actual surface 
profile measured by AFM. The last curve presents the 
maximal value of H c % that was calculated at a given angle 
for profiles measured by the AFM. In case B, when the 
surface is relatively smooth H c3 /H C 2 = 1.56 for the paral- 
lel to the surface dc field. This is in good agreement with 
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FIG. 8: (Color online) . The angular dependence of H c s/H C 2 
ratio: experiment and calculations for several models of the 
surface roughness. The curve A corresponds to the ideally 
smooth surface; the curve B - to sinusoidal surface with m = 
20 and n = 1; curve C - to the model surface with m — 10 
and n = 5; and curve D - to the measured by AFM profile. 
See text. 



the experimental value 1.57. However, the calculated H c3 
decreases with an angle and for 6 — ir/2 H c3 /H C 2 = 1 
while the experimental value is 1.28 (Fig. [8]). For a more 
steep roughness H c3 /H C 2 ~ 1.25 at angles close to 7r/2 
(curve C). This ratio fits well the experimental value of 
1.28 for a perpendicular dc field. However, the curve C 
was calculated for a very steep roughness that was not 
observed for our sample by AFM. We have to note that 
the scale at this model rough surface is several tenths of 
nanometers. One could expect the comparatively large 
H c 3 in the normal dc field if there are sufficient number 
of spots at the surface for which the angle between the 



dc field and the local normal is considerably large, but 
the AFM data do not confirm this. We see that actual 
surface roughness could not yield such a large value of 
i? C 3 / H C 2 in a perpendicular dc field. Another hypothesis 
that there is a superconducting network in the bulk which 
is not seen by dc measurements [3] at Ho > H C 2 leads to 
the requirement that the bulk has to have the normal 
conductivity that is approximately 100 times larger then 
was measured. 



V. CONCLUSIONS 

We have studied the low frequency ac response of a 
YB 6 single crystal in tilted dc magnetic fields. The de- 
veloped approach has permitted us to obtain the value of 
the ac surface current at each side of the sample. The H c3 
angular dependence when dc field rotates out of the plane 
was obtained. The experiment showed that in the per- 
pendicular to the surface dc field the transition to super- 
conducting state took place in a field Hq « 1.28i7 C 2. This 
unexpected experimental result cannot be ascribed to 
the possible affect of the surface roughness. We demon- 
strated also that the induced ac surface current is a func- 
tion of an angle between ac and dc fields when both fields 
are parallel to the surface. 
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